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(57)Abstract: 

PROBLEM TO BE SOLVED: To estimate the intake air flow introduced in a 
combustion chamber through an intake passage precisely regardless of an 
engine operation state. 

SOLUTION: The calculation device 30 of an internal combustion engine 10 
estimates the temperature after the prescribed term from the present 
point of time of the intake air staying on the down streamside part 28 than 
the throttle valve 24 of the intake passage 20 based on the fluid model 
regarding the intake air in an intake passage 20 and the temperature 
information at the present point of time of the intake air. The calculation 
device 30 estimates the intake air flow introduced in the combustion 
chamber 29 from the downstream part 28 after the prescribed term from 
the present point of time while reflecting the estimation result. 
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CLAIMS 



[Claim(s)] 

[Claim 1] In the inhalation air-content prediction equipment of the internal combustion engine which adjusts the amount 
of the inhalation air introduced into a combustion chamber by the throttle valve of an inhalation-of-air path To the 
inhalation air in the aforementioned inhalation-of-air path The related fluid model The temperature information on this 
time of this inhalation air The amount of the inhalation air introduced into the aforementioned combustion chamber 
from the aforementioned downstream portion in a predetermined period from this time, presuming the temperature after 
a predetermined period from this time of the inhalation air which is alike, is based and exists in the portion of a 
downstream rather than the aforementioned throttle valve of the aforementioned inhalation-of-air path, and reflecting 
this presumed result Inhalation air-content prediction equipment of the internal combustion engine characterized by 
having a prediction means to predict. 

[Claim 2] In the inhalation air-content prediction equipment of an internal combustion engine according to claim 1 the 
aforementioned prediction means The fluid model built based on the heat energy conservation law and mass 
conservation law about inhalation air of the aforementioned downstream portion, Based on the pressure and temperature 
in this time of inhalation air of the amount in this time and temperature of the inhalation air which flows into the 
aforementioned downstream portion, and the aforementioned downstream portion, the pressure and temperature after a 
predetermined period are presumed from this time of the inhalation air of the aforementioned downstream portion. 
Inhalation air-content prediction equipment of the internal combustion engine characterized by being what predicts the 
aforementioned inhalation air content based on the this pressure presumed and temperature. 

[Claim 3] In the inhalation air-content prediction equipment of an internal combustion engine according to claim 2 the 
aforementioned prediction means mt (0) and temperature for the amount in this time of the inhalation air which flows 
into the aforementioned downstream portion Tha, Pm (0) and temperature for the pressure in this time of the inhalation 
air of the aforementioned downstream portion Tm (0), When the ratio of specific beat of mc (i) and air is set to kappa 
and capacity of R and the aforementioned downstream portion is set [ the amount of the inhalation air introduced into a 
combustion chamber from the aforementioned downstream portion ] to V for a gas constant, The pressure Pm (i) of the 
inhalation air of the discrete formula (1) obtained from each aforementioned conservation law, (2) and the 
aforementioned inhalation air content mc (i), and the aforementioned downstream portion, and relational expression 
with temperature Tm (i) (3) 
[Equation 1] 

P m W = rp mV ~ ») + ■ * ■ y ' K ( ' " l) ' T ha " m c<* ' l > ' T m« - »)) "(l) 
-p2-(0 = ^-0 - 1) + A/ . £ . («,(/ - i> - m c (t - i)) • • -(2) 
m c<» = f{F m (t),T m V)) —(3) 

At.Mftmm t/wd = £(a«j)) 

Inhalation air-content prediction equipment of the internal combustion engine characterized by being what presumes the 
pressure Pm (n) and temperature Tm (n) after the predetermined period tfwd from this time of the inhalation air of the 
aforementioned downstream portion by performing processing only whose number of times n of predetermined it is 
alike and repeats the based operation. 
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[Claim 4] In the inhalation air-content prediction equipment of an internal combustion engine according to claim 3 the 
aforementioned prediction means Perform a series of aforementioned repetitive operation for every predetermined 
interruption period, and it sets to this interruption timing. Inhalation air-content prediction equipment of the internal 
combustion engine characterized by being what computes either [ at least ] the pressure Pm in this time of the inhalation 
air of the aforementioned downstream portion used for the repetitive operation in next interruption timing (0), or the 
temperature Tm (0). 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[The technical field to which invention belongs] This invention relates to the inhalation air-content prediction equipment 
of the internal combustion engine which predicts the amount of the inhalation air introduced into an engine combustion 
chamber. 
[0002] 

[Description of the Prior Art] In case the amount (inhalation air content) of the air introduced into the combustion 
chamber of an internal combustion engine through an inhalation-of-air path determines the controlled variable of 
internal combustion engines, such as fuel-injection control and ignition-timing control, as what reflects engine 
operational status most with an engine rotational frequency, it is used, and this is usually called for based on the 
detecting signal from an air flow meter or an intake-pressure sensor prepared in the inhalation-of-air path. 
[0003] By the way, it faces computing the controlled variable of an internal combustion engine based on an inhalation 
air content, and if it is at the time of transient operation when the inhalation air content in this time called for from 
detecting signals, such as an air flow meter, is used as it is, this controlled variable comes to be determined based on a 
different inhalation air content from the time of the controlled variable actually being reflected in engine control. Thus, 
if the inhalation air content used for the determination of a controlled variable differs from the inhalation air content in 
case the controlled variable is reflected in engine control, even though control adapted to actual engine operational 
status cannot be performed but it will raise the control precision, a limitation comes to arise naturally. 
[0004] Then, the inhalation air content after a predetermined period is predicted from this time, and the equipment 
which determined the controlled variable based on the forecast is conventionally proposed so that it may be indicated by 
JP,2-42160,A. By performing such prediction, the inhalation air content at the time of a controlled variable being 
reflected in engine control can be foreseen beforehand, this controlled variable can be determined now, and 
improvement in control precision comes to be achieved. 

[0005] Moreover, an inhalation air content changes also with an engine rotational frequency, or intake-air temperatures 
besides the opening of a throttle valve. For example, since the inhalation density of air becomes large so that an intake- 
air temperature becomes low, an inhalation air content also comes to increase. Therefore, it becomes important, when 
obtaining a prediction result accurate [ also taking into consideration the influence of such an intake-air temperature ], in 
case an inhalation air content is predicted as mentioned above. 

[0006] The intake-air temperature in the position of an upstream is detected, and it is made to make the intake-air 
temperature (for this to be hereafter, called "upstream intake-air temperature 11 , when it is necessary to distinguish 
especially from the intake-air temperature in other positions of an inhalation-of-air path) reflect in prediction of an 
inhalation air content rather than the throttle valve of an inhalation-of-air path with this point and the above-mentioned 
equipment. 
[0007] 

[Problem(s) to be Solved by the Invention] However, although the above-mentioned conventional equipment is making 
the intake-air temperature reflect in prediction of an inhalation air content, it cannot but perform the prediction to the 
last by making for position-change and its time change of an intake-air temperature not to arise during the prediction in 
an inhalation-of-air path into a prerequisite. 

[0008] If surely it is at the time of steady operation by which the opening of a throttle valve is held uniformly, since it is 
in the range which most position-change of the intake-air temperature in such an inhalation-of-air path does not have, 
and can also disregard change of a time intake-air temperature, though an inhalation air content is predicted only based 
on the upstream intake-air temperature mentioned above though the inhalation air content was predicted under the 
above-mentioned prerequisite, gross errors do not arise in the prediction result. 
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[0009] However, if it is when the opening of a throttle valve changes rapidly at the time of transient operation, intake-air 
temperatures come to differ according to the position in an inhalation-of-air path, and the time change cannot be 
disregarded greatly, either. Therefore, the above-mentioned prerequisite does not meet the actual temperature state of 
the inhalation air in such an inhalation-of-air path. 

[0010] For example, when the amount of the inhalation air which flows into the downstream portion of this valve when 
the opening of a throttle valve increases increases rapidly, the inhalation air of this downstream portion will be 
compressed and the temperature (henceforth a "downstream intake-air temperature") will come to rise more temporarily 
than the upstream intake-air temperature mentioned above. Since it will be conversely expanded by the inhalation air of 
this downstream portion when the amount of the inhalation air which flows into the aforementioned downstream portion 
on the other hand when the opening of a throttle valve decreases decreases rapidly, a downstream intake-air temperature 
comes to fall more temporarily than an upstream intake-air temperature. 

[001 1] thus, temporary, if it is at the time of transient operation from which the opening of a throttle valve changes - an 
imitation - a downstream intake-air temperature may become temperature which is changed sharply in time and is 
different from an upstream intake-air temperature Consequently, if it is in the conventional equipment which predicts by 
making for change [ be / time / position-/ it ] of an intake-air temperature not to arise in an inhalation-of-air path into a 
prerequisite, so to speak, the influence by the dynamic behavior of the intake-air temperature at the time of such 
transient operation is not reflected in prediction of an inhalation air content, therefore the error of the prediction result 
cannot be disregarded, either. 

[0012] This invention is made in view of such the actual condition, and the purpose is in offering the inhalation air- 
content prediction equipment of the internal combustion engine which can predict the amount of the inhalation air 
introduced into a combustion chamber through an inhalation-of-air path with a sufficient precision even if it is at the 
transient operation time. 
[0013] 

[Means for Solving the Problem] The means and its operation effect for attaining the above-mentioned purpose are 
indicated below. In the inhalation air-content prediction equipment of the internal combustion engine which adjusts the 
amount of the inhalation air with which invention indicated to the claim 1 is introduced into a combustion chamber by 
the throttle valve of an inhalation-of-air path To the inhalation air in the aforementioned inhalation-of-air path The 
related fluid model The temperature information on this time of this inhalation air The amount of the inhalation air 
introduced into the aforementioned combustion chamber from the aforementioned downstream portion in a 
predetermined period from this time, presuming the temperature after a predetermined period from this time of the 
inhalation air which is alike, is based and exists in the portion of a downstream rather than the aforementioned throttle 
valve of the aforementioned inhalation-of-air path, and reflecting this presumed result It is made to have a prediction 
means to predict. 

[0014] A prediction means presumes the temperature of the inhalation air of the aforementioned downstream portion 
based on the fluid model about the inhalation air in an inhalation-of-air path, and the temperature information (physical 
quantity which has these or correlation, such as a mode of temperature or a temperature change) on this time of this 
inhalation air, and makes the presumed result reflect in prediction of an inhalation air content in the above-mentioned 
composition. Therefore, even if it was the case where the temperature of this downstream portion changed from this 
time transitionally, when the transitional change has been grasped, an inhalation air content can be predicted, and a very 
accurate prediction result can be obtained now. 

[0015] Moreover, since the amount of the inhalation air introduced into a combustion chamber through an inhalation-of- 
air path tends to increase, so that the pressure of the inhalation air of the above-mentioned downstream portion becomes 
high, and so that the temperature becomes low, it can ask for this based on the above-mentioned pressure and 
temperature. Therefore, the amount of the inhalation air introduced into a combustion chamber in a predetermined 
period from this time based on the pressure and temperature which are they-presumed can be predicted now by 
presuming the pressure and temperature after a predetermined period from this time of the inhalation air of a 
downstream portion. 

[0016] Furthermore, the pressure and temperature after a predetermined period can presume these based on the heat 
energy conservation law about the inhalation air of a downstream portion, and a mass conservation law by relating with 
physical quantity called the pressure and temperature of inhalation air in this time of the amount in this time and 
temperature of the inhalation air which flows into a downstream portion, and its downstream portion from this time of 
the inhalation air of this downstream portion. 

[0017] Invention indicated to the claim 1 is set to the above-mentioned inhalation air-content prediction equipment as 
composition materialized more like invention which indicated to the claim 2. therefore, the - aforementioned prediction 
means The fluid model built based on the heat energy conservation law and mass conservation law about inhalation air 
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of the aforementioned downstream portion, Based on the pressure and temperature in this time of inhalation air of the 
amount in this time and temperature of the inhalation air which flows into the aforementioned downstream portion, and 
the aforementioned downstream portion, the pressure and temperature after a predetermined period are presumed from 
this time of the inhalation air of the aforementioned downstream portion. The composition that it is what predicts the 
aforementioned inhalation air content based on the this pressure presumed and temperature is employable. 
[0018] Invention indicated to the claim 3 is set to the inhalation air-content prediction equipment of an internal 
combustion engine according to claim 2. moreover, the aforementioned prediction means mt (0) and temperature for the 
amount in this time of the inhalation air which flows into the aforementioned downstream portion Tha, Pm (0) and 
temperature for the pressure in this time of the inhalation air of the aforementioned downstream portion Tm (0), When 
the ratio of specific beat of mc (i) and air is set to kappa and capacity of R and the aforementioned downstream portion 
is set [ the amount of the inhalation air introduced into a combustion chamber from the aforementioned downstream 
portion ] to V for a gas constant, The pressure Pm (i) of the inhalation air of the discrete formula (1) obtained from each 
aforementioned conservation law, (2) and the aforementioned inhalation air content mc (i), and the aforementioned 
downstream portion, and relational expression with temperature Tm (i) (3) 
[0019] 
[Equation 2] 

P m P> = P mV - J > + A ' * * * y ' { m t <' " ] > • T ha " m c^ - 0" I^V - ' *0) 
P P R 

■^r-d) = -=r" (* - I) + Af • — ' { m t (' - 0 ~ ™ C V - d) ' ' '(2) 
*c<0 = /(^(0,7>)) ---(3) 

f{P m «>, T m a)): P m ii), T m wfco i, vc© ■» 

By performing processing only whose number of times n of predetermined it is alike and repeats the based operation, it 
is supposed that it is what presumes the pressure Pm (n) and temperature Tm (n) after the predetermined period tfwd 
from this time of the inhalation air of the aforementioned downstream portion. 

[0020] Since according to the above-mentioned composition it is dispersed and expressed as the aforementioned fluid 
model is shown in each above-mentioned formula (1) and (2), the pressure Pm (n) and temperature Tm (n) after a 
predetermined period can be presumed based on a comparatively easy algorithm from this time of the inhalation air of 
the aforementioned downstream portion. 

[0021] It sets to the inhalation air-content prediction equipment of an internal combustion engine according to claim 3 
like invention indicated to the claim 4 in performing such data processing, moreover, the - aforementioned prediction 
means Perform a series of aforementioned repetitive operation for every predetermined interruption period, and it sets to 
this interruption timing. If the composition that it is what computes either [ at least ] the pressure Pm in this time of the 
inhalation air of the aforementioned downstream portion used for the repetitive operation in next interruption timing (0) 
or the temperature Tm (0) is adopted Not using a sensor etc., ** can also ask for the pressure Pm of the above- 
mentioned inhalation air (0), or its temperature Tm (0) serially through the above-mentioned repetitive operation, and 
simplification of the composition of inhalation air-content prediction equipment can be attained now. 
[0022] Moreover, generally, since the responsibility of a general-purpose temperature sensor is low, when the 
temperature Tm of inhalation air (0) is detected using this sensor, it is difficult to detect change of this temperature Tm 
(0) with a sufficient precision, and it cannot finish being avoided that an error arises in the detection result. 
[0023] In this point and the above-mentioned composition, if especially the temperature Tm of an inhalation air content 
(0) is computed, while being able to attain simplification of composition, the fall of the predictability resulting from the 
detection error of such a sensor can be avoided. 
[0024] 

[Embodiments of the Invention] The 1st operation form of this invention is explained below [the 1st operation form]. 
[0025] It is made to apply the inhalation air-content prediction equipment concerning this invention to the internal 
combustion engine equipped with an air flow meter as a sensor which detects an inhalation air content with this 
operation form. First, with reference to drawing 1 , the outline of this internal combustion engine 10 and the 
composition of this prediction equipment are explained. 

[0026] The inhalation-of-air path 20 of an internal combustion engine 10 is constituted by the suction-port 23 grade 



1 



Page 4 of 14 



which connects the surge tank 22 connected to the aisleway of the throttle body 21 in which the throttle valve 24 was 
formed, and the downstream of this throttle body 21 and this surge tank 22, and a combustion chamber 29. After 
metering of the inhalation air which flows the inhalation-of-air path 20 is carried out by the throttle valve 24, it is 
introduced in a combustion chamber 29 at the time of valve opening of an inlet valve 26. 

[0027] Inhalation air-content prediction equipment predicts the amount of the inhalation air introduced in a combustion 
chamber 29 in a predetermined period from this time, and is greatly constituted by the various sensors 41-45 including 
an arithmetic unit 30 and an air flow meter 41 . 

[0028] The arithmetic unit 30 is equipped with the operation part 31 which performs data processing, and the memory 
32 the various function data used on the occasion of the program concerning such data processing or its execution were 
remembered to be. This arithmetic unit 30 is constituted by the electronic control which performs various control of an 
internal combustion engine 10. 

[0029] Others, an intake temperature sensor 42, the throttle sensor 43, the cam angle sensor 44, and crank angle sensor 
45 grade are contained in sensors 41-45. [ air flow meter / above-mentioned / 41 ] All, the air flow meter 41 and the 
intake temperature sensor 42 are formed in the portion of an upstream from the throttle valve 24 at the inhalation-of-air 
path 20, and detect the amount and temperature of the inhalation air which flows in a surge tank 22 through this valve 
24 (henceforth "the upstream intake-air temperature Tha"), respectively. 

[0030] Moreover, with this operation form, the thing of the heat ray formula equipped with the detecting element 410 
which consists of a heat ray as the above-mentioned air flow meter 41 is used. Drawing 2 shows the cross-section 
structure of this detecting element 410. As shown in this drawing, this detecting element 410 is constituted by the heat 
ray 412 which consists of platinum etc., and the glass layer 414 which covers the circumference of this heat ray 412. 
[0031] The throttle sensor 43 is formed near the throttle valve 24, and detects the opening (throttle opening TA) of this 
valve 24. The crank angle sensor 45 is formed [ near / where the cam angle sensor 44 carries out the opening-and- 
closing drive of the inlet valve 26 / the cam shaft (illustration abbreviation) ], respectively near the crankshaft 
(illustration abbreviation) which carries out the rotation drive of this cam shaft. An arithmetic unit 30 incorporates the 
output signal of these cam angle sensor 44 and the crank angle sensor 45, and detects the rotational speed (engine 
rotational frequency NE) of a crankshaft, its rotation phase (crank angle CA), and the opening-and-closing time (valve 
timing VT) of an inlet valve 26, respectively. In addition, the mechanism (illustration abbreviation) in which the 
opening-and-closing time of the above-mentioned inlet valve 26 is changed based on engine operational status, such as 
the engine rotational frequency NE and an engine load, is prepared in the internal combustion engine 10 of this 
operation form. Therefore, the above-mentioned valve timing VT will be suitably changed according to engine 
operational status. 

[0032] Next, prediction processing of the inhalation air content by the equipment of such this operation form is 
explained. The principle which will be the requisite for this prediction processing first is explained. 
[0033] The space from the inhalation-of-air path 20 to a combustion chamber 29 from a throttle valve 24 First, the 
portion of an upstream From 27 and this throttle valve 24, rather than an inlet valve 26 by the downstream The portion 
of an upstream (It is hereafter called an "upstream portion") (It is hereafter called a "downstream portion") It divides 
into 28 and three portions called a combustion chamber 29, and the fluid model about the inhalation air which exists in 
the downstream portion 28 among each [ these ] portions 27-29 is built. 

[0034] Namely, the amount of the inhalation air which flows out of "M" and this downstream portion 28 per time the 
mass of the inhalation air which exists in this downstream portion 28 in a combustion chamber 29 If the amount 
(henceforth a "throttle- valve through put") of the inhalation air which flows (it is hereafter called "the inhalation air 
content in a cylinder") into the downstream portion 28 from "mc" and the upstream portion 27 is set to "mt" Variation 
[ per time of the above-mentioned mass M ] d (M)/dt can be expressed like the following formula (5) based on a mass 
conservation law. 
[0035] 
[Equation 3] 

-jf M = m t - m c •"•(5) 

On the other hand, the heat energy, in which this inhalation air has them when "Tm", its isochore specific heat, and the 
isotonic specific heat are set to "Cv" and "Cp", respectively serves as "Cv-M-Tm" in the temperature (henceforth a 
"downstream intake-air temperature") of the inhalation air which exists in the downstream portion 28. Moreover, since 
the heat energy which the inhalation air which flows into a combustion chamber 29 from "Cp-mt-Tha" and the 
downstream portion 28 has [ the heat energy which the inhalation air which flows into the downstream portion 28 from 
the upstream portion 27 has ] is "Cp-mc-Tm", the variation d (Cv-M-Tm)/dt per time of above-mentioned heat energy 
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Cv-M-Tm can be expressed like the following formula (6) based on a heat energy conservation law. 

[0036] 

[Equation 4] 

j;{C v - M T^-Cp mr Tha-Cp mc T* " '(6) 

Moreover, about the capacity of the downstream portion 28, if the equation of state about the inhalation air which exists 

in the downstream portion 28 sets the pressure (henceforth a "downstream intake pressure") of the inhalation air for "V" 

and the said division 28 to "Pm" and sets a gas constant to "R", it will become like the following formula (7). 

[0037] 

[Equation 5] 
P m V=M-A-T n —(7) 

Here, if the ratio of specific beat of air is set with "kappa" (=Cp/Cv), a formula (5) can express the following formulas 
(9) and a formula (6) like the following formulas (8) using the above-mentioned formula (7). 
[0038] 
[Equation 6] 

= «■ y-imfTho-mc Tm) ••■(*) 

■ f •*■•--) ■•'"» 

Furthermore, if a discrete interval is set to "deltat", the left part of each [ these ] formula (8) and (9) can be expressed as 
shown in the following formula ( 1 0) and ( 1 1 ). 
[0039] 
[Equation 7] 

in - Pm (O ~ PmV- '> ... fl0) 

Therefore, the following discrete formulas (12) and (13) are drawn from each [ these ] formula (8) - (1 1). In addition, in 
this formula (12) and (13), "i" in a parenthesis shows a value and "i-1" shows a value last time this time, respectively. 
[0040] 
[Equation 8] 

= P m V-l) + At-K-y- («,</ - i) - T ha - m c o - i) • r m( / - i))- • -(12) 
P P R 

Y~ (/) = O - i) + Af . — . {m t (i-n- m c o -})) • • -(13) 

On the other hand, the inhalation air content mc in a cylinder can be calculated using the following empirical formula 
(14). 
[0041] 
[Equation 9] 
m c = ISL { a .p m -b) •• # (14) 

In an upper formula (14), each of "a" and "b" is constants, and is set as the flow of the inhalation air in the portion from 
the downstream portion 28 to a combustion chamber 29 at the affecting engine control parameter, i.e., here, based on the 
engine rotational frequency NE and valve timing VT. The relation between these constants a or a constant b and the 
engine rotational frequency NE, and valve timing VT is beforehand defined by experiment, and is memorized by the 
memory 32 of an arithmetic unit 30 as a function map. Moreover, "To" is the reference temperature (constant) of the 
inhalation air when setting up each constants a and b by such experiment. 

[0042] As shown in each following formula (15) and (16), if this empirical formula (14) is used, each above-mentioned 

formula (12) and (13) can deform. 

[0043] 

[Equation 10] 
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P»m = P n V - i) + A/ • K y ■ (m,u - l) • r te - T 0 ■ (a ■ P„« - D - *)) ■ ■ -(1 5) 
^-(0 = ^0 - i) + A/-£.[m,<; - ,) - 7 JSL_. (a ■ P„ V - i) -ft)] -(16) 

Moreover, it can ask for the downstream intake-air temperature Tm using the formula (17) of the following obtained 

from each [ these ] formula (15) and (16). 

[0044] 

[Equation 11] 
T m (» = -^- -07) 

771 



Each [ these ] formula (15) So that clearly from - (17) The throttle- valve through put mt in this time, It is based on the 
downstream intake pressure Pm, the downstream intake-air temperature Tm, and the upstream intake-air temperature 
Tha. The downstream intake pressure Pm and the downstream intake-air temperature Tm after the above-mentioned 
discrete interval deltat can be presumed after a predetermined period (i.e., here) from this time. Furthermore, by 
substituting for a formula (14) the downstream intake pressure Pm presumed in this way and the downstream intake-air 
temperature Tm, the inhalation air content mc in a cylinder after discrete interval deltat can be presumed from this time. 
[0045] Moreover, it is made to calculate the above-mentioned throttle-valve through put mt in prediction processing of 
this operation form according to two ways the followings differ. Hereafter, the method at the time of calculating this 
throttle- valve through put mt is explained. 

[0046] First, the procedure at the time of calculating the throttle- valve through put mt based on the detecting signal of an 
air flow meter 41 is explained. As everyone knows, an air flow meter 41 detects the amount of said using the heating 
value taken from this heat ray 412 changing according to an inhalation air content, in case inhalation air passes the 
circumference of a heat ray 412. The relation between the detecting signal of this air flow meter 41 and an inhalation air 
content is beforehand called for by experiment, and is memorized by the memory 32 of an arithmetic unit 30 as a 
function map. Therefore, an inhalation air content can be calculated based on the detecting signal and this function map 
of an air flow meter 41 . 

[0047] However, as compared with a heat ray 412, generally, since heat capacity is large, the above-mentioned glass 
layer 414 has the inclination for the temperature Tgs to be overdue to change of the amount of the inhalation air which 
passes near the detecting element 410, and to change. For this reason, if there is an inhalation air content calculated from 
the above-mentioned function map if it is at the time of steady operation from which an inhalation air content does not 
change at the time of transient operation from which a true inhalation air content and this inhalation air content of a 
match change, the inhalation air content calculated from a function map comes to shift from a true inhalation air content. 

[0048] Then, in order to foresee beforehand the response delay about the temperature change of such a glass layer 414 
and to calculate a true inhalation air content, the thermal equilibrium model about this layer 414 is built, and it is made 
to compute an inhalation air content with this operation gestalt based on the reverse model. In addition, in the following 
explanation, a true inhalation air content is made into the throttle-valve through put mt, and the inhalation air content 
acquired based on the detecting signal and the above-mentioned function map of an air flow meter 41 is distinguished as 
a map operation value GA. 

[0049] First, if temperature of a heat ray 412 is set to "Tht" (refer to drawing 2 ), the variation d (Tgs)/dt of the 
temperature Tgs of the glass layer 414 can be expressed like the following formulas (18) from the balance with the 
heating value transmitted to the glass layer 414 from a heat ray 412, and the heating value transmitted to inhalation air 
from this layer 414. 
[0050] 

[Equation 12] 

A ' Tt T & = B ( Tht " T *) " ( C + D ' ' l T " " Th °) "* (18) 

Here, "A", "B", "C", and "D" are constants determined according to **, such as the cross section of a heat ray 412, 
length and its resistivity, and a heat transfer rate between the glass layer 414 and a heat ray 412, a heat transfer rate 
between the glass layer 414 and inhalation air. 

[0051] Moreover, at the time of steady operation, since transfer of the heat between the glass layer 414, and a heat ray 
412 and inhalation air is lost, as shown in the following formula (19) and (20), the left part d (Tgs)/dt of temperature 
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changes of a formula (18), i.e., the amounts of the glass layer 414, is set to "0", and the map operation value GA and the 

throttle- valve through put mt become equal. 

[0052] 

[Equation 13] 
GA = mt ~<20) 

Therefore, the following formula (21) is obtained from each above-mentioned formula (19) and (20). 
[0053] 

[Equation 14] 



Furthermore, if a formula (18) is transformed using this formula (21), the following formula (22) from which the 

temperature Tgs of the glass layer 414 was eliminated will be obtained. 

[0054] 

[Equation 15] 



Moreover, the left part of this formula (22) can be expressed like the following formula (23), if a discrete interval is set 

to "deltat." 

[0055] 

[Equatio n 16] 

d4GA _ yjGAtj) - jGAu - i> 
dt to "" l } 

Therefore, the following discrete formulas (24) are drawn from each above-mentioned formula (22) and (23). In 
addition, in this formula (24), "alpha" and "beta" are constants which become settled by the above-mentioned constant 
A-D, and "j" in a parenthesis shows the value and "j-1" shows the value last time this time, respectively. 
[0056] 

[Equation 17] 



The throttle-valve through put mt in this time (j) can be presumed based on the map operation value GA in this time (j), 
and the map operation value GA in front of discrete interval deltat so that clearly from this formula (24). 
[0057] Next, the procedure at the time of calculating the throttle-valve through put mt based on engine operational 
status, such as the throttle opening TA, is explained. When a part for opening of a throttle valve 24 is assumed to be a 
kind of orifice, the amount (throttle- valve through put mt) of the fluid which passes this orifice is expressed by each 
following formula (25) and (26) as everyone knows. 
[0058] 

[Equation 181 






■*(P m ,P ff ) -.-(25) 




Here, "mu" is the flow coefficient of an orifice (a part for opening of a throttle valve 24), and "Pa" is atmospheric 
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pressure (constant value). Moreover, "At" is the effective-area product of a throttle valve 24, and can ask for this 
uniquely from the throttle opening TA. 

[0059] Here, if the time of steady operation by which the throttle opening TA is held uniformly is assumed, a formula 

(25) will become like the following formula (27). 

[0060] 

[Equation 19] 

«iw ■ M -At ■ -7=%= • •fc.w. p a) "-(27) 

Here, "mtTA" and "PmTA" show the throttle-valve through put mt at the time of steady operation, and the downstream 
intake pressure Pm, respectively. The throttle- valve through put mtTA at the time of these steady operations and the 
downstream intake pressure PmTA can be determined as the flow of the inhalation air in the portion from the upstream 
portion 27 to a combustion chamber 29 at the affecting engine control parameter, i.e., here, based on the throttle opening 
TA, the engine rotational frequency NE, and valve timing VT. The relation between the throttle- valve through put mtTA 
at the time of these steady operations or the downstream intake pressure PmTA, the throttle opening TA and the engine 
rotational frequency NE, and valve timing VT is beforehand called for by experiment, and is memorized by the memory 
32 of an arithmetic unit 30 as a function map. 

[0061] If this formula (27) is used, a formula (25) can be expressed like the following formula (28). 
[0062] 

[Equation 20] 

Mt - m,TA (28> 

Furthermore, if it is at the time of steady operation, since the throttle- valve through put mtTA and the inhalation air 
content mc in a cylinder become equal (mtTA=mc) and the downstream intake-air temperature Tm and the upstream 
intake-air temperature Tha become equal (Tm=Tha), the above-mentioned formula (14) can be expressed like the 
following formula (29). 
[0063] 

[Equation 21] 

<n<TA = * ?mTA " *) " '<2*> 

'ha 

And if this formula (29) is used, the above-mentioned formula (28) can be expressed still like the following formula 

(30). 

[0064] 

[Equation 22] 

- f£r-(- ■ - •)• - (30) 

The throttle- valve through put mt can be presumed based on the value PmTA of the downstream intake pressure Pm 
when assuming that it is at the upstream intake-air temperature [ in this time ] Tha and downstream intake-pressure Pm, 
atmospheric pressure Pa, and steady operation time (TA, NE, VT) so that clearly from this formula (30). 
[0065] It is made to predict the inhalation air content mc in a cylinder after a predetermined period to be the throttle- 
valve through put mt calculated based on each option in this way from this time with this operation form based on each 
formula (14) - (17) mentioned above. In addition, in the following explanation, the throttle-valve through put mt asked 
for the throttle-valve through put mt calculated based on engine operational status, such as the throttle opening TA, 
based on the detecting signal of the throttle-valve through-put estimate mtl and an air flow meter 41 is distinguished as 
throttle- valve through-put measured value mt2. 

[0066] Next, drawin g 3 and drawing 4 are combined, referred to and explained about the detail of the prediction 
processing by the equipment of this operation form. Drawing 3 and drawing 4 are flow charts which show an example 
of the procedure at the time of predicting the inhalation air content mc in a cylinder. This processing of a series of is 
performed with an arithmetic unit 30 the interruption period (every [ For example, " ] 8msec. M ) for every predetermined 
time. 

[0067] On the occasion of this processing, the prediction period tfwd at the time of predicting an inhalation air content is 
set up first (Step 100). This prediction period tfwd will be the fuel oil consumption in the last fuel injection, and a thing 
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set up according to fuel injection duration, if it puts in another way. Therefore, the inhalation air content mc in a 
cylinder after the time which the last fuel injection took from this time passes will be predicted here. Incidentally, this 
fuel injection duration is computed through processing concerning fuel injection with this another processing. 
Moreover, in processing concerning this fuel injection, fuel oil consumption is computed based on the newest inhalation 
air content mc in a cylinder computed through this processing. 

[0068] next, the downstream intake pressure Pm after [ this time to ] the prediction period tfwd, the downstream intake- 
air temperatures Tm, and those ratios -- Pm/Tm is computed based on the aforementioned throttle- valve through-put 
estimate mtl (Step 200) In addition, the downstream intake pressure Pm after the prediction period tfwd computed 
below based on the throttle-valve through-put estimate mtl the downstream intake-air temperatures Tm and those ratios 
- Pm/Tm with "Pmltfwd", "Tmltfwd", and "Pml tfwd/Tmltfwd", respectively moreover, the downstream intake 
pressure Pm in the middle of the calculation, the downstream intake-air temperatures Tm, and those ratios - it writes 
"Pml (i)*\ "Tml (i)", and "Pml/Tml (i)", respectively [ Pm/Tm ] 

[0069] On the occasion of calculation of this downstream intake-pressure Pmltfwd and downstream intake-air 
temperature Tmltfwd, processing only whose number of times of predetermined repeats the operation based on each 
above-mentioned formula (15) - (17) and a formula (30) is performed. 

[0070] That is, as shown in drawin g 4 , the increment of the counter value i showing the number of times of a loop of 
this processing (initial value : "0") is carried out first (Step 202), and it is judged whether next the prediction periods 
tfwd are deltatl or more predetermined times (Step 204). Here, the above-mentioned predetermined time deltatl is set 
up equally (for example, "8msec. 11 ) to the interruption period of this routine. And when the prediction period tfwd is 
judged to be these deltatl or more predetermined times (Step 204: YES), discrete interval deltat of each above- 
mentioned formula (15) and (16) is set up equally to a predetermined time deltatl (Step 206). 
[0071] next, based on a formula (30), the throttle-valve through-put estimate mtl (i) computes - having - further - 
after that and each formula (15) - (17) - being based - the downstream intake pressure Pm 1 (i), the downstream 
intake-air temperatures Tml (i), and those ratios - Pml (i)/Tml (i) is computed, respectively (Step 208) And a 
predetermined time deltatl is subtracted from the present prediction period tfwd, and the subtraction value (tfwd-delta 
tl) is set up as a new prediction period tfwd (Step 210). 

[0072] Next, when it is judged whether this new prediction period tfwd is below "0" (Step 212) and this prediction 
period tfwd is judged to be longer than "0" (Step 212: NO), processing after Step 202 is performed again. Under the 
present circumstances, if it is judged that the prediction period tfwd is less than deltatl predetermined time in Step 204 
(Step 204: NO), discrete interval deltat will be set up equally to the prediction period tfwd (Step 207), and processing 
after Step 208 will be performed. 

[0073] And by repeating processing of Steps 202-212 If the prediction period tfwd becomes below "0" (Step 212: YES) 
the downstream intake pressure Pm 1 (i), the downstream intake-air temperature Tml (i), and a ratio -- Pml/Tml (i) of 
them downstream intake-pressure Pmltfwd after the prediction period tfcvd, downstream intake-air temperature 
Tmltfwd, and a ratio - it is set up as Pml tfwd/Tmltfwd, respectively (Step 214) 

[0074] For example, when the interruption period of this routine is "8msec." and the prediction period tfwd is set as 
"30msec." through processing of previous Step 100, "Times [ 4 ]" Discrete interval deltat will be set up in order through 
processing of the above-mentioned step 204,206,207 in the mode "8msec." ->"8msec." ->"8msec." -> "6msec", and 
processing to Step 202 - Step 212 will be repeated, and the downstream intake pressure Pm 1 (4) in case the counter 
value i is "4", the downstream intake-air temperatures Tml (4), and those ratios Pml/Tml (4) ~ downstream intake- 
pressure Pmltfwd after [ this time to ] the prediction period tfwd, downstream intake-air temperature Tmltfwd, and a 
ratio ~ it will be set up as Pml tfwd/Tmltfwd, respectively 

[0075] Then, while the counter value i was reset by "0", after discrete interval deltat is set up equally to the above- 
mentioned predetermined time deltatl (Step 216), processing is returned to Step 300 shown in drawing 3 . 
[0076] At Step 300, the map operation value GA in this time (j) and (it is hereafter written as "GA(0)") are computed 
based on the detecting signal of an air flow meter 41 . Next, based on the map operation value GA in this this time (0), 
and the map operation value GA (j-1) (it is hereafter written as "GA (-1)") computed in the last interruption timing, the 
throttle- valve through-put measured value mt2 in this time (0) is computed using a previous formula (24) (Step 400). 
[0077] and this throttle- valve through-put measured value mt2 (0) - being based - each previous formula (15) - (17) - 
using » the downstream intake pressure Pm after discrete interval deltat (i), the downstream intake-air temperatures Tm 
(i), and those ratios - Pm/Tm (i) is computed, respectively (Step 500) In addition, the downstream intake pressure Pm 
computed in this step 500 (i) the downstream intake-air temperatures Tm (i) and those ratios - the value computed to 
this interruption timing among Pm/Tm (i) with "Pm2 (1)", "Tm2 (1)", and "Pm2/Tm2 (1)", respectively It writes "Pm2 
(0)", "Tm2 (0)", and "Pm2/Tm2 (0)", respectively. [ the value computed to the last interruption timing ] 
[0078] each following formula (31), (32), and (33) - being based - the downstream intake pressure Pmtfwd after [ this 
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time to ] the prediction period tfwd, the downstream intake-air temperatures Tmtfwd, and those ratios - while 
Pmtfwd/Tmtfwd is computed, the inhalation air content mctfwd in a cylinder after the prediction period tfwd is 
computed by being based on the formula (34) obtained from the above-mentioned formula (14) (Step 600) [ next, ] 
[0079] 

[Equation 23] 

Pmtfrd = ?m\tfwd - + P*2«» " # ( 31 ) 

Pjntfwd ^mXifvfd P m \ P m o 

~ — = -=^(0) + (0) •••(32) 

l mtfwd l m\tfwd 'mi l ml 

^--5*^ -(33) 
r m tjwd 

Tmtfwd 
F mtfwd 



In this way, after the inhalation air content mctfwd in a cylinder is computed, the initial value of data processing in next 
interruption timing is set up, respectively (Step 700). 

[0080] Drawing 5 is a timing chart which shows the setting mode of such initial value. In this drawing, while engine 
operation is started in time tl, this routine shall be started. 

[0081] Initial value [ in / the interruption timing (time t2 and t3) of the 2nd henceforth / as shown in this drawing ], 
Namely, the downstream intake pressures [ Pm / Pm (0) and / 2 ] 1 (0), the downstream intake-air temperature Tml (0), 
Tm(s)2 (0), these ratios Pml/Tml (0), Pm2/Tm2 (0), and the map operation value GA (-1) It is set up as follows using 
the value computed in the last interruption timing as the values (Pml (1) etc.) or the value (GA (0)) in this time after 
discrete interval deltat, respectively. 
[0082] 

[Equation 24] 

Pml (0) +- Pml (I) 
Pm2«» «- PmlV) 
T m \(0) <" T m \0) 

^(0)^-^0) 
GA{-\) <- GA(Q) 

On the other hand, about the initial value in the first interruption timing (time tl), it is set up as follows from the 

calculation value in the last interruption timing not existing. 

[0083] 

[Equation 25] 

Pml CO) «- Pa 
Pm2(°) «- Pa 
Tmim «- Tha 
Tml «0 <" Tha 

P m\ Pa 

l m\ i ha 

- p- 

GA{-\) <r~ 0 



Thus, after setting up each initial value in next interruption timing, processing of this routine is once ended. 

[0084] As explained above, with this operation form, a fluid model is built based on the heat energy conservation law 

and mass conservation law about inhalation air of the downstream portion 28. The fluid model, Based on the throttle- 
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valve through put mt, the upstream intake-air temperature Tha, the downstream intake pressure Pm, and the downstream 
intake-air temperature Tm, the downstream intake pressure Pmtfwd after the prediction period tfwd and the downstream 
intake-air temperature Tmtfwd are presumed from this time. Furthermore, it is made to predict the inhalation air content 
mctfwd in a cylinder after the prediction period tfwd from this time based on them. 

[0085] (1) Therefore, even if it was the case where the downstream intake-air temperature Tm changed from this time 
transitionally, when the transitional change has been grasped, the inhalation air content mc in a cylinder can be 
predicted, and a very accurate prediction result can be obtained now. 

[0086] In case the fluid model of such inhalation air is built, the mass conservation law and heat energy conservation 
law about inhalation air of the downstream portion 28 (2) Moreover, a previous formula (12), It is made to presume the 
downstream intake pressure Pmtfwd after the prediction period tfwd, and the downstream intake-air temperature 
Tmtfwd from this time by a discrete formula's as shown in (13) expressing, and performing processing which repeats 
the operation based on these dispersion formula (12) and (13) the number of predetermined times. Therefore, the 
downstream intake pressure Pmtfwd required for prediction of the inhalation air content mctfwd in these cylinders and 
the downstream intake-air temperature Tmtfwd can be presumed now based on a comparatively easy algorithm. 
[0087] (3) In case the above-mentioned repetitive operation is performed a predetermined interruption period, it is made 
to compute the downstream intake pressure Pm (0) and the downstream intake-air temperature Tm (0) in this time used 
to next interruption timing in this interruption timing furthermore. For this reason, not using the sensor for detecting the 
pressure and temperature in the downstream portion 28 etc., ** can also ask for these downstreams intake pressure Pm 
(0) and the downstream intake-air temperature Tm (0) serially through the above-mentioned repetitive operation, and 
can attain now simplification of the composition of inhalation air-content prediction equipment. 
[0088] (4) By presuming the downstream intake-air temperature Tm in this time (0) especially, the need of generally 
using a general-purpose temperature sensor with low responsibility is lost, and the fall of the predictability resulting 
from the detection error of such a sensor can be avoided now. 

[0089] Moreover, it is made to calculate the throttle-valve through-put estimate mtl presumed from engine operational 
status, such as the throttle opening TA, and the throttle- valve through-put measured value mt2 computed based on the 
detecting signal of an air flow meter 41 as a throttle-valve through put mt with this operation gestalt, respectively. 
[0090] Here, when the downstream intake pressures [ Pm / Pm (i), Pm 12 I (0), and / 2 ] 1 (1) called for using these 
throttle- valves through-put estimate mtl and the throttle-valve through-put measured value mt2 and the downstream 
intake-air temperature Tml (i), Tm2 (0), and Tm2 (1) compare with an actual intake pressure and an actual intake-air 
temperature, it is checked by experiment of this invention persons that there are the following inclinations. 
[0091] That is, about the downstream intake pressure Pm 1 (i) called for using the throttle- valve through-put estimate 
mtl, the downstream intake-air temperatures Tml (i), and those ratios Pml/Tml (i), there is an inclination for the time 
change to reflect correctly an actual time change of an intake pressure or an intake-air temperature especially. On the 
other hand, about the downstream intake pressures [ Pm / Pm (0) and / 2 ] 2 (1) called for using the throttle- valve 
through-put measured value mt2, the downstream intake-air temperature Tm2 (0), Tm(s)2 (1) and those ratios Pm2/Tm2 
(0), and Pm2/Tm2 (1), there is an inclination for it to be well in agreement with an intake pressure, an intake-air 
temperature, or the actual value of those ratios at the time of steady operation. 

[0092] Fortius reason, if it is at the time of transient operation by computing the downstream intake pressure Pmtfwd 
after the final prediction period tfwd, and the downstream intake-air temperature Tmtfwd using each previous formula 
(31) - (33), while being able to make a time change of an actual intake pressure and an intake-air temperature reflect in 
the calculation correctly, it is at the time of steady operation, [0093]. [ i.e., ] 

P equation 26] 
mtXfwd = PmiW 

P mt\fwd _ P m \ 
miXfwd 1 m\ 

It is [0094] when ******(ing). 
[Equation 27] 

Pmtfwd = Pml (0) 
Pmtfwd P m 7 

— = -=™-<o) 

s mtfwd l ml 



From a bird clapper, the downstream intake pressure Pmtfwd and the downstream intake-air temperature Tmtfwd come 
to be called for only based on the throttle-valve through-put measured value mt2. 
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[0095] (5) therefore - according to this operation gestalt - the time of steady operation - be - the time of transient 
operation - be ~ the downstream intake pressure Pmtfwd after the prediction period tfwd and the downstream intake-air 
temperature Tmtfwd can be correctly presumed from this time, as a result the inhalation air content mctfwd in a cylinder 
can be predicted more now to accuracy 

[0096] (6) Furthermore, with this operation gestalt, since the thermal equilibrium model about the glass layer 414 of an 
air flow meter 41 is built and it is made to calculate the above-mentioned throttle- valve through-put measured value mt2 
based on the reverse model, the response delay about the temperature change of this layer 414 can be foreseen 
beforehand, and it can ask for this. Therefore, the more exact throttle-valve through-put measured value mt2 can be 
calculated, and the downstream intake pressure Pmtfwd and the downstream intake-air temperature Tmtfwd can be 
presumed now still more correctly based on this. 

[0097] The 2nd operation gestalt of this invention is explained focusing on difference with the 1st operation gestalt 
below [the 2nd operation gestalt]. 

[0098] With this operation gestalt, the point it is made to apply the inhalation air-content prediction equipment 
concerning this invention to the internal combustion engine equipped with an intake-pressure sensor as a sensor which 
detects an inhalation air content is different from the 1st operation gestalt. As shown in drawing 6 , this intake-pressure 
sensor 46 is attached in the surge tank 22, and detects the pressure Pm of the inhalation air of the interior, i.e., a 
downstream intake pressure. 

[0099] Next, the detail of the prediction processing by the equipment of this operation gestalt is explained. Drawing 7 is 
a flow chart which shows an example of the procedure at the time of predicting the inhalation air content mc in a 
cylinder. This processing of a series of is performed with an arithmetic unit 30 the interruption period (every [ For 
example, ,f ] 8msec") for every predetermined time. In addition, at Step 100 shown in this drawing 7 , and Step 200, 
since the same processing as Step 100 and Steps 200-216 which are shown in previous drawing 3 and previous drawing 
4 is performed, explanation of the content of processing is omitted. 

[0100] After each processing of these steps 100,200 is performed, at Step 350, the downstream intake pressure Pm in 
this time (i) and (it is hereafter written as "Pm3(0) t! ) are computed based on the detecting signal of the intake-pressure 
sensor 46. Next, the downstream intake pressure Pm detected in the downstream intake pressure Pm 3 (0) detected in 
this interruption timing, and the last interruption timing (i-1) (it is hereafter written as "Pm3 (-1)") Based on the 
following formula (35) obtained from a previous formula (15), the throttle-valve through put (to namely, last 
interruption timing) mt (i-1) (henceforth "the throttle- valve through-put measured value mt3 (-1)") in front of discrete 
interval deltat is computed from this time (Step 450). 
[0101] 

[Equation 28] 

„,,«_,> = _l_ j^o)-/W-> .JL. + To ( a . M - »)} .-OS) 



Further And this throttle- valve through-put measured value mt3 (-1) and a previous formula (16), It is based on the 
following formula (36) obtained from (17), and (37). The downstream intake-air temperature Tm in this time (0) The 
ratios Pm3/Tm3 (0) of (it being hereafter written as ,! Tm3(0)"), and the downstream intake pressure Pm 3 (0) and the 
downstream intake-air temperature Tm3 (0) are computed, respectively (Step 550). 
[0102] 

[Equation 29] 
£» 



C0) = Al<-» + A, . £ . ^ /3( -i) - • (a - />„ 3 <-i) - A)] - -(36) 



^-4^ ^(37) 

7 m3 



In addition, in the above-mentioned step 450 and Step 550, when computing the downstream intake-air temperature 
Tm3 (0) etc. based on a formula (35) and (36), the value which the value into which the upstream intake-air temperature 
Tha was read to the last interruption timing was used, and was set up based on the last engine rotational frequency NE 
and last valve timing VT in interruption timing about each constants a and b is used. Moreover, when this time is the 
first interruption timing, the value in this time [ intake-air temperature / upstream / Tha ] and the value to which each 
constants a and b are set based on the engine rotational frequency NE (= "0") and valve timing VT at the time of an 
engine halt are used. 

[0103] Next, the formula which replaced "Pm2 (0)" and "Pm2/Tm2 (0)" by "Pm3 (0)" and "Pm3/Tm3 (0)" in each 
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previous formula (31) and (32), respectively, each formula (33) and (34) - using -- the downstream intake pressure 
Pmtfwd after [ this time to ] the prediction period tfwd, the downstream intake-air temperatures Tmtfwd, and those 
ratios, while Pmtfwd/TmtfWd is computed The inhalation air content mctfwd in a cylinder after the prediction period 
tfwd is computed (Step 650). 

[0104] In this way, after the inhalation air content mctfwd in a cylinder is computed, the initial value of data processing 
in next interruption timing is set up, respectively (Step 750). Drawing 8 is a timing chart which shows the setting mode 
of such initial value. In this drawing, while engine operation is started in time tl, this routine shall be started. 
[0105] Initial value [ in / the interruption timing (time t2 and t3) of the 2nd henceforth / as shown in this drawing ], The 
downstream intake pressures [ Pm / Pm (0) and / 3 (-1) ] 1, the downstream intake-air temperature Tml (0), Tm(s)3 (-1) 
and these ratios Pml/Tml (0), and Pm3/Tm3 (-1) [ namely, ] It is set up as follows using the values (Pm3 (0), Tm3 (0), 
etc.) calculated in the last interruption timing as the values (Pml (1) etc.) or the value in this time after discrete interval 
deltat, respectively. 
[0106] 

[Equation 30] 
r m i<o) - r w id) 

7^3 <-»><- 7^3 «» 
tn\ 'ml 

Pm3 M><- P M ^m 

On the other hand, about the initial value in the first interruption timing (time tl), it is set up as follows from the last 

calculation value or last detection value in interruption timing not existing. 

[0107] 

[Equation 31] 
K\W <- Pa 

Pm2i-»<- Pa 
T m \ (0) - Tha 

P P 
Pm*i-V<- Pa 

Thus, after setting up each initial value in next interruption timing, processing of this routine is once ended. 
[0108] Also in the 2nd operation gestalt explained above, the operation effect of an abbreviation EQC can be acquired 
about (1) indicated in the 1st operation gestalt, (2), and (4). moreover, about the operation effect of (3) indicated in the 
1st operation gestalt, conversely [ of that for which the intake-pressure sensor 46 for detecting the downstream intake 
pressure Pm is needed separately ], since an air flow meter 41 does not need the sensor which becomes unnecessary and 
detects the downstream intake-air temperature Tm at least, at this point, simplification of the composition of inhalation 
air-content prediction equipment can be attained - it becomes like 

[0109] Furthermore, the downstream intake pressures [ Pm / Pm (0) and / 3 ] 3 (1) called for based on the detecting 
signal of the intake-pressure sensor 46, the downstream intake-air temperature Tm3 (0), Tm(s)3 (1) and those ratios 
Pm3/Tm3 (0), and Pm3/Tm3 (1) have the inclination for it to be well in agreement with an intake pressure, an intake-air 
temperature, or the actual value of those ratios at the time of steady operation. Therefore, the operation effect equivalent 
to (5) indicated in the 1st operation gestalt also according to this operation gestalt can be done so. 
[0110] Each operation gestalt explained above can also change and carry out composition as follows. 
- Although the estimate based on engine operational status, such as the throttle opening TA, and the actual measurement 
based on the detecting signal of an air flow meter 41 or the sensor of intake-pressure sensor 46 grade are calculated as 
an amount of the inhalation air which passes a throttle valve 24, respectively and the inhalation air content mc in a 
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cylinder was predicted with each above-mentioned operation gestalt using each [ these ] value, based on either these 
estimate or an actual measurement, you may be made to perform this prediction. 

[01 1 1] - Although it was made to predict the inhalation air content mc in a cylinder with each above-mentioned 
operation gestalt supposing the internal combustion engine 10 equipped with the mechanism in which the valve timing 
of an inlet valve 26 is changed, also in the internal combustion engine which is not equipped with this mechanism, the 
inhalation air content mc in a cylinder can be predicted according to the same procedure. 

[01 12] - When the mechanism which controls the flow of the inhalation air of the inhalation-of-air path 20, for example, 
the swirl valve which controls the intensity of the swirl generated in a combustion chamber 29, performs the above- 
mentioned prediction again in the internal combustion engine formed in this inhalation-of-air path 20, set up the above- 
mentioned constants a and b, the downstream intake pressure PmTA at the time of steady operation, and the throttle- 
valve through put mtTA as a function of the controlled variables (for example, opening of a swirl valve etc.) of the 
above-mentioned mechanism. 

[01 13] - In prediction processing of each above-mentioned operation gestalt, although the empirical formula (14) which 
makes this air content mc the function of the downstream intake-air temperature Tm and the downstream intake pressure 
Pm was used when calculating the inhalation air content mc in a cylinder, the method of calculating this air content mc 
is not limited to this. 

[01 14] - Although the downstream intake pressure Pm and the downstream intake-air temperature Tm were computed 
with each above-mentioned operation gestalt by dispersing each above-mentioned formula (8) obtained as a fluid model 
of inhalation air, and (9), and performing a loop operation, you may make it compute these downstreams intake pressure 
Pm and the downstream intake-air temperature Tm based on the operation expression which solved each [ these ] 
formula (8) and (9) analytically. 

[01 15] - Although the above-mentioned prediction period tfwd was set as the value suitable for calculation of this fuel 
oil consumption with each above-mentioned operation gestalt that the inhalation air content mc in a cylinder used for 
calculation of fuel oil consumption should be predicted, this prediction processing is applicable also to other engine 
control, such as ignition-timing control, by setting up the prediction period tfwd suitably. 

[01 16] - Although the interruption period of a prediction manipulation routine shown in drawing 3 or drawing 7 and 
discrete interval deltat were equally set up with each above-mentioned operation gestalt, this discrete interval deltat can 
also be changed "1/2 of range which can determine initial value in next interruption timing, i.e., interruption period," 
twice etc., for example so that it may become 1/n time (n : integer). 



[Translation done.] 



Page 1 of 5 



* NOTICES * 

Japan Patent Office is not responsible for any 
damag s caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



DRAWINGS 




[Drawing 1] 



31 



10 -MM 29-tm 

20 -MM 30-HRt 

24-20 9 V i* 41-1750^-9 

28-T1M* 42-mrj* 



[Drawing 2] 



410 




414 412 



[Drawing 3 ] 



Page 2 of 5 



( » 


r } 






ifflBffcfwdt* 



BftflBtfwdftOTWBttB'u. 



Zast-MUMTsTNHSA (0) «f B 



-S200 



■^S300 



7»7S»»GA (-1) . SA (0) CMUT 

Mftwoji-jWfiamftite (o) eta 



ttra&tto-Rtt&ip* (i) 

. TfiUBST* (1) . J&WT* (i) sw 



•S4O0 



-S500 



P* (0) *-P« (1) 
P-<0) ♦FW(i) 
T* (0) (1) 
T-(0)*T*(1) 
P-/T u (0) «P*/T* (1) 
P«/T^(0)*P-i/T*(l) 
GA (-1) <-GA (0) 



-S600 



-S700 



c 



3 



[Drawing 5] 



IMMttol 



mm 



p«(i) ■» 




P*(i) i 




p-(l) 


p«(l> 




FU(1) 




p«(i) 


T^(i) 




T-(l) 




T„(l) 


T*(l) 






T«(l) 






T*(lT 


lu 






^(1) 

lei 






In 














§»(D 
la 


GA(0) - 






GA(0) - 




GA(0) 


P.4 (0) 




'P*(0) 




r P«ttt» 


P«(0) 




P«(0) 




P«(0) 


T*(0) 




T„«)) 




Tift) 


T*(0) 




T*(0) 


-> • 


T*(0) 










ft" 










I ft" 


GA (-1) 




GA(-i) 1 


u 


GA (-1) 



ti 



t2 



t3 



it 



At 



[Drawin g 4] 



( y -S200 



-S202 




YESJ06 



At*4ti 



S207 



At«-tfwd 



mu (i) . ft» (i) . T„ (i) . ^ (i) ©M f" S208 



T 



tfwd*tfwd-Atl -S210 




tfwd<0 

YES 



S212 
NO 



P*tt*a*Pm. (i) 
T«««d*Ti (i) 

P^tfMl/TBItfa^Pu/TM (i) 



■-S214 



I 



i«-rOJ H-S216 
At«-Ati 



( w-o ) 
[Drawing 6] 



21 2B 20 

2 7 24 I ri 1 ' ^44 



on 




[Drawing 8] 



Tu (1) 



P«(0) 
T«(0) I 

*« n 



p. 



L (0) 



P«(0) 
T«(0) 

ft 

P-(-l) 



ti 



P-(l) 
T-(i) 



P»(0) 
T«(0) 
P. 



To 



(0) 



r P- (0) 
T*«» 

^(0) 

ltd 

Po(-l) 
T*(-l) 



P-(i) 
T«(l) 



.(0) 
.(0) 



ft" 

P« (0) 
T«(0) 

• u 

rP-(-l) 
T«(-l) 

I S3 



t2 



t3 









r 4t * 


At 





[Drawing^] 



) 



1 




nmtfwdK -sioo 









-S200 



tMft«(-l) f P-<0)CWl»I _ 
tUfflAtHMDy hA^GWKiikB (-1) «W 



WiTOTHflSfflT- (0) . 

J&WU(0)«l 



WBMtfwdta 
Rn«ASUm<t M C«iii 



P M (0)*P.,(1) 
Ti(0)«-T..(l) 
P-/T* (0) ♦P^/T- 11) 
P*(-1)*P«(0) 
T«(-1)*T«(0) 
P«/T*(-1)«-P«/T*(0) 



-S450 



--S550 



-S650 



--S750 



V 



C & ) 



Page 5 of 5 



[Translation done.] 



f&m 2001-41095 

(P 2001-41095 A) 

(43)4^0 ¥J5£13^2^13 0 (2001.2. 13) 



(51) Int. CI. 7 SSSUE-JI- 
F 0 2 D 45/00 3 6 6 

41/18 



F I 

F 0 2 D 45/00 
. 41/18 



f^-r (##) 
3 6 6 Z 3G084 
F 3G301 



»3E«* *»* W*K*>«4 OL 



51) 



<2i)a«## 


#^¥11-215535 


(71) HUB A 


000003207 














(22)ttl«B 


1*£7J1 29 0 (1999. 7. 29) 












(72)3§PJ?# 


/h# A^ 










**D«*ffl rfi h 3 ^BTTl#it& 


h3^gib 
















(72)*W# 


Ag W 












h 3^g®, 
















(74) ABBA 


100068755 




















1 





(54) imw<D*m AKtinosASXifist 



(57) [KM] 

[MHl BR«3i»S:ai:r««fe^rt«c»AStLS®A^ 




/ 10-flMW 29-ttti 

31 20-ftSl* 30 

24-2D»hJl* 41-l7^D^-9 



(2) 



«FBB 2001-41095 



[flt*« 2 ] fl|#£ 1 E*Ort«»HOKA3S«*^a!l 
flc^/ui:, 89ETffiffly«P^«-«fcA-4-S»A^«om«F 



KA^ftoS^^^SfSMBIHS^ffi^Xr^ia**:*^ 
U K«t**lxSBEAatftaaE«J:*<5v^TmE»Affi«; 

Affi»*^H3g«. 
[W*5C3] »*«2«cE«ort««H<ogRAffiiS[*^ 

WE^JH^att. «ftETaE«»4>«waKAi-5»A2SftO 
aW-iSl?«>**:mt(0)XtfaL«*rTha. HuETSEffltt^ 

(o). METK«A^&«Htt(c«ASix5ftAS«; 
^t^mc(i), ffiftGJJfcfRJtSric. m*^i5:^R. ME 
TaB*B»4^©S««rV t Ufci:#. jN*E#ft#Iija>&ft 
p>H««»a; (D , (2) . 3feWcii9ERA£SC*m 
c(i) fc^TaE«as»<0!RAS«Offi*Pni(i)at/ia* 
Tm(i)fcOBB«* (3) 
[SCI] 



-■>•**■ -/n c a -i> r m o-i)>--(i) 

= ±3L U _ 1)+ A/ (m t o- \)-m c « - i)) • • -(2) 
tn nt 

m c (i) = f(P m iiUT m V)) ---O) 



/(^« 0), T m <#)): (#>, r OT (o [r o v *x<7> B3& 



(flw/ = £ (A* «) 



C^cJ;^ «reT«««ro»«5»Affl«oa«P^bBf 30 
SSBPfltfwdlfcCDJE^ Pm(n)^OTSTm(n) $r«^i"5 fc 
otfcs: £ t««i-r5rtJB«BB<oiKAffi«*^«K 

Mo 

[»**4l W**3|dE«^rt«S«H«>KAS**^ 

S9ET^ftflffl^KA^«om^'ST^JE*Pm(0)Rt^ 
?aSTffl(0)©^< t ^-^$:StUi-^^^T'fe^^i: 40 
4r»« <b r 5 rtil^H©!ft A£« J»«««. 

[00 0 1] 

WA*n5»A^«o*^^j|8r-5rt«S«H<OKAffi« 
[0 0 0 2] 

A^tt^^moft (»AS«f) f4^«BBlHWBRti:t> 



[0 0 0 3] tr5T% «Affi«*KS^Trt«&«B8 

«§B*J»^*»lcR!fcSixS«FjftLi:tt»ft5»Affl** 

(^J^iS^^I-ffi^/cKA^fi^:, ^<o*J»i*s«t 
BaffiSP l-R^: £ ti Z> t ft cdpJ£ A3£«ft £ 4ot^5 

[0 0 0 4] ^^T% ^fl?2-4 2 1 6 0*4*«I-E 
tt£*i5J:5K:, ««Fjft^e>9f^n«<OBRA35«** 

o^tX\ 0Jff*^«H0Jff{cSK$^SR«>KAS« 



3 

[0 0 0 5] »A3£«*tt* «BBlHlte«[^ay 

tf , K«S.^ffi < -5 15 < 4 5 

[0 0 0 6] r©jS, ±SB3£«rti\ RftiB^Dy 

5) Sr»Affi**^WcR»S-frSJ:5«c:urv>S. 
[0 0 0 7] 

*<^ag«tt. p^A^^o^ay{ce^^u^^^^T^ 

br^<0^«l£:?T5 tor 
[0 0 0 8] BfcWc, * n y h/l'flWBBAas— JEWS!* 

t sm-c # 5 ffiH tb z> it *b s ±.mmm^w 

fc±^flJK«a^^£^#KA^ft4^^JWS: 
[0 0 0 9] fcC^3ftS, i§atWEI*lC*5V^T^i3y 

^o^^atc^b-rs^^feorfi. wtfaimm 

lEttftfefrlii 9 U^!»SCiiKrtic*5itS»Affi«;o 

[ooio] , * n v b/is&cDfflmfrmttz z 

t\c£*j, l^#(0T«L«jffl5>^«EA-rs»AS«O*35S 
bfc»§\ I^TSKflUffl^OKA^^jEESIS 

KiJ: 19 . fifrlBT«IE«?B^^i«Ai"SKASao» s A» 
KSM> Ufc^lcteig^ ^T«E«l«»o»AS«^Bg 

[0 0 11] rod: 9{w^n 5/ h/u#oBH«36S^l:-f 5 
««»6B*K:*>ortt:* — H^ttfcHirJ:. Tfl£««®SU&# 

V'&Zo *<o&m. »«iiKrt^*5u>r»«aoffiBW 

$rfT9^*co^g^feor^ r 5 vtzmmmumc^s 



(3) 2001-41095 

4 

[0012] z-<D5£W&^ o b/cHtf^^r^^^yc 
twfeo, ^o@WttK«SBS»i§i:r*Rfig^rt«-* 
A£*t*»A^flw>**. ifi«jite«Frfeorfc»ffifi* 

< ^SJ-TS rt^H0»A£9l&rFfl!J£K 
[0 0 13] 

OKASStft^SJKfiJcifcv^r, iMEKfiuiKrt^KA 
^i-^i"^^^^/u^fi]p^A^(7)^^r^^S 
W«i:^S^^rmifaK«iiK(OffiriB^ci y 9 

^5»A^«0*Sr^ai-t"S^ai#aSr«XLSJ; 5ia 
20 rv^>5 0 

[0 0 14] ±SB*i£K*5^r. ^m^&te. K«5S8S 
rtOKA^«tcBlli-Si«»^/^tl^®A^OS^S- 

KA^^SrlSL, ^*3e*g*&KA£«*o 
^SIKRSfcSitrSo Sot, ^coTii»(DK^I 

ftSEftSriGa u^±r® Asaasr^si-r s r t fcX 
Mftxm&<D&\,^mM$k&'&z>zt&x%z>£? 

30 \Z/j:^ a 

[0 0 15] ^fc. K«iiKS:iiDr*S«feSli:*AStt 
5»Affi«tO*i, -hgSTMJSP^OKA^om^^ 

ifts^fcSfcJ?). rft£±SEJ£fc;RT*aSKS^^r#to 
5:i#«5. fot, T*ffi(ffi#<&«A^*03BH# 

f?. ^tL^^^nSEE^^^{Cg<5^r^.^^ 
^^fflra«JcJ»«|^fc»A$HSKA2ESloSSr^ai 
i"5 r t fcx # s J: 9 (c/i^ 0 

40 [0 0 16] r©TMiJS5)^®A^(DSM 

»OKAS«(-ov^ro^R^/^=¥^R#lUMW^flfl:^S 
#9J^-S<5v^r. T«LiBja5»lc:iffiAi-SKA3g«toa^F 

[0 0 17] ttor, f9*«ltcfE<SU^:*l8S:J:<9* 
MUcK^LT, lf*3S2lcfe®Ufc^0^(DJ: o 
\^ ±EKAffi«*^SJSSI»c*3V^r, -SWE^W^S 
50 S(rlEr«EfflJffl»OKA^«(cB8i-5f»^^^^r«# 



(4) 



2001-41095 



met. ^^mmTmm\n^(D^A^%(D^^(DEE 
[0018] *tz. m^m3\cmm^t^m^ 



* ^T^t$rmt (0) RXfUm. £: Tha. mZTVrM&ftcD 
» A^U<D^^ X <Dl£t) &Pm(0)R tfffi.m $rTm(0), 
«MET«ffil»#a> e>«**^*A $ tt* « A^tDML* 
mc(i). 35«01t«Ut*: jc , «ffc£S:«:R, fl9i2TfflEffl!J 

(1) , (2) . afeWdBuEKA^SCftmcCpi: 
miST«EfiiJffl^oKA^ftWJE^Pm'(i)Xt52afiTm(i) 
ir^K^ (3) 
[0019] 
10 [|&2] 



=P «C- »> + • * • — . (/*,(/ - l) - 7^ - ot c (/ - 1) • T m {i - l))- • •(!) 



CO = If- U - i) + Af ~ • (m, </ 
trt 



i)-m c (i-\)) "-(2) 



OT c<0 = f( P mOhT m (n) ---(3) 



/ ( p m T m <o): P m V), T m v *T o Rft 



«JlET8ltfl!Jffl^^»A^ft;om«F^d»e>9f 
^fflMtfwdfg^JE^ Pni(n)M^STni(n) £#£5£-r5 
<OX&>Z>b LT^£ 0 

[0020] ±fiB«^:fcj:tLtf. WiaaEfr^e^/udsjblE 
(l) , (2) {^*tL5J:5{wie»^$tLT*S 
^T^5lt^ HulET«Kfti]«a5»oeRASSLom«FJ#C*^ 

^^>o 30 

[0021] r. 5 LfcSW&3&fT5K*>fc*K 
««Bi^»A^«*^aiJ3S«^i3v^-c, -ffirE^sj^a 

fJii^-f S v^^<^»S«L3S^^:fflv^fettSfl^rl^T^KfllJSP 
#CD® A^S^TOfOEA Pm(0)M£STjn(0) CO 

Sffli-tLf*. ±IE®A^SKDJE*Pm(0). Sc^f«:tf>t& 
aTmCOSr-fe^^ESrfflv^ttiE^fcjg&SSriiCT 40 
®Affi«ft^»JS»^«fi8<z>fB 

[0 0 2 2] £fc. infflcoiaffi-t^-^n-t^ 

«4a*ffi^fc«>* Pir^SrfflV^KA^SlOiaSTni 
(0)Sr«tH-r5<t 5(-L/cS^ Hi&«Tm(0)<7>»ffcSr 

[0 0 2 3] ±ia*jSfcio1/*-(\ 4*K»Affi* 

i(DSgTin(0) SrHCttJ-T 5 «fc 5 JC-rtltf . ffijSL<Dffimt 
]Z>ZkfcX%Z>b kh\^. Z 9 Ltz±>-V-<DWiltim 50 



H^gEg-r S^aSfiJSofiTSrlEliB'r 5 r 1 &X% 5 J: 
[0 0 2 4] 

com l o*lBeilSJcov^TlftW"t'S. 
[0 0 2 5] KA^a**«ltBi-S-fe 

\^x^T7u—*—?*ffi^%iftffimm\^&w 

[0 0 2 6] F^*?S«ii 1 0 cORfuISg 20(t * * S> h 
/u#2 4#Rltbhfc^Py h/UsKv^ 2 1 cop^aJil 
rco^n-/ h/Ustf^-f 2 1 (OT^flJtC^RSttfclJ- 

££&Jfc^*RSt#— h 2 3^KJ:9«^£tLT^5o 
»£ui&2 OSrStttSRA^SCti:, ^oyWu#2 4|: 
J:oTWftStL*:8L K«#2 6<OBB#^«S8fe^2 9 
rt^WASttS. 

[oo2 7] KA^^s^aijKfifL a^^&BfSffl 

Bfl»tt:*&#5S2 9^tc^A*tLSPAA^cofi$:^9J-r 

4 1 $rt^Di6i:i-5^S-tr>'-th4 l~45^[cJ;o 
T«^$ttTt^o 

[oo2 8]-»jiRt3oii, mm*&m&nn-rz>Wi& 
tftc[gLTffiv>fetLS#aii^--^^fets$ttfc^^- 

Xffif&£ixZ> 0 
[0 0 2 9] ir>^~4 1—4 5^fi. ±.m^TZ7a — * 



7 

^tl^o i77B-^-?4 1, "M^yf4 2li 

^Tttfc, ft«UiB2 O ic:^v>T^Dy h/u#2 4 J: t) 
±8Kffllo8P»{cKJtP>tbTio9. HW2 4 ^btt- 

[0 0 3 0] £fc, *HifeJglli-ett. ±|Ex77P-^ 

^4U lt. «y»^&*s«iasfP4 i o &ffi*_itm 
&&.<Di > <D&m^btix^z>* B2ii C<D$l/±iS5 4 1 

mas 4 lod 0»jx.*ie^«E^p>*s!Ry»4 i2t, ^ 
^4 1 2<Dmm*ffi&tz>tf7*m4 i 4 taot 

[0 0 3 1] y h/U-fe:>i^4 3lt y h/U#2 

4<Z>iS«^RJtbHT*5 9. I^#2 4^K (*n-/h 
«TA) $:^ttJi-^o ^7A^ir>1^4 4ti\ 
2 6S:HEH^»i-5*A^^7 h (BI^BS) <offi«fc:. 

bt^TV^5o ^$1^3 Ofi. rttb^7^^-ir>1h4 4 
^^^V^^-^^4 5CD{±J;^fi#£&9i£^ *7 
^7 hOl§Mtei$ffi («BB!eME*NE) RX**<om 

mm (^7>^ca) »i:pm#2 6op»i 

SStt^ffi<Drt«*«BI 1 0 tCH _LlEPj*ja# 2 6 (DMEflB* 

±ffi^^^-f ^ v^vTt>«B85itettffiiciEi:r 
[0032] r o Lttmmmm<Dmmc£.z>»& 

[0 0 3 3] £i\ 8RftiiK2 0 d^JjRSgS 2 9^X'0 

r±8Eifi!jgP^J £V^5) 27, 13* n h/U-#2 4 £ «9 
T8E«J-CK«# 2 6 X V i>±9M<o«6& C£AT> rTifft 
fllffll^J £^5) 2 8. MSi2 9^^f:3oO 
S5#£#«\ rtL^>#a55>2 7-2 9 60 9 t.TMiJM 

[0 0 3 4] EP*>, rOT8Kfi!ja5»2 8icff&-f5KA 
^Ofi$r I'm j . ^TSEfflUa55>2 8^e>«Btttt2 9 

»A£«*J £V^50 £ r mc j . ±SEffiiJ«P5J-2 7^5> 
TSE<WSB» 2 8 ic:8(EA-t-SKA^«^ft («T\ n 
y £ r m tj t-f%t\ _hfE® 

SMOH#PdiSfc «9 £>g£ffcft d (M) / d t ttKftft^FBJ 
lc£^T&5£ (5) W5(cSt^^^5o 



(5) ftBH 2001-41095 

8 ' 

[0 0 3 5] 
[13] 

d 

— A/ = mt - m c 

TSK«Sffl5»2 8lC#«E-rSKA^m^i&S (J£* 

mRXJ^EEitmtr^tt^ti revj , repj b-tZk* 
|SJ®A^«^^T"t-Sf»^^/U^ftt TCv • M • Tmj t# 
So £fc, Jb^{»Jffl5>2 7^e>T^««»5>2 8J^iKA-t- 
10 SP^A^C0^"i-Sl^^/^^ rep • mt • Thaj , 

•f-^Sia^/n^flS rCp-mc-TmJ X'frZ^kfrb. 

(CvM'Tm) /dt tt*fc^*/l'*«#lU»-S-3l^T 

(6) (0«J:9lcttw<l:^x^5 0 . 
[0 0 3 6] 

[14] 

™(Cw - ^ • r m ) = C P - mt Tha- C p mc T m 

20 T«M««#2 8{^ffii-SKASSC«cov^O* 

ft*S*J4. TJSfEfiiJ6iJ3)'2 8 0«aS: Tvj , m&&2 
8X*<Di&7<&n<0&J3 (J^T. fT«fEfiiJ»ftU=Ej 

9) ^ rp m j N trj 

(7) (OX 5l-*5o 
[0 0 3 7] 
[15] 

p m -V = M ■ R- Tm ~>(7) 

::t\ &%i<Dttmtt$: r Kj (=c P /c v ) tio< 

JilH* (7) £r/BV>T5S (5) fil^TO^ (9) x 
30 * (6) liJ^ToS; (8) oj:5^tri:^T^6 0 
[0 0 3 8] 
[16] 

Ztlh^ (8) , (9) tfXfci2te> gttfcr^Pi^ 
fa t j b-fZ>h&L#: do), (ii) t^j: 

40 [0 0 3 9] 
[17] 

.^■1 ff'"-7f"-" ..,,„ 
Air. J A' c 1 

^or. ^ttb^^: (8) ~ (ii) ^P>KTO«l»S; 
(12), (13) ^«#UJ$ixSo ft. -«o^: (i 
2) , (13) tdioVNT, r i j tt^EMfi. 

50 f i - l j fiHulHlffiSr^ix^ix^i- 0 



[0 0 4 0] 



(6). 



* * [f£8] 



mm 2001-41095 



10 



Oo = ^(/-i) + Arv |-(ffi f (,-i).7i fl -« c (/-i)-r w (#- i)>--(i2) 

(13) 



P P R 

Y~ to = ^-0 - i) + A/ ~ • (/n, <f - i) - m c (> - i)) 



[0 04 1 ] 
[&9] 



ft^llKS: (14) 



- fe) ---(14) 



_t5£ (14) tci^r. r a j . 

T*£> 0 . T$L®mft 2 8 d^&Mttt£ 2 9 £T-<OgJ5#tc::fc 



fi^ftH^tci: 9j£fct>iv SfM30W^ey 3 2 
iCiift^yytLTlEtt^n-CV^So £fc, l"Toj 

10 «A£«<DS2pia* (J£%0 -CfcS. 

[0 0 4 2] :(D^ (14) Srffll^Si:. _bfE#5£ 
(12), (13) It &<£>#i£ (15), (16) CO 

[0 0 4 3] 
[f&l 0] 



-co = ^-o - i) + to 



R ( 

y ^W f (i-l)- 



.<'-»-*)) 



•(16) 



TSEffl»£li&Tm«: % ^thh^ (1 5) , (1 
6) ^P>»P>*t5«Toa (17) S:l^t*Wri: 

[0 0 4 4] 
[ftl 1] 



T m iO = 



0> 



ci s) — (i 7) m 

TaE<B»«fflTm*r-a; (14) fcftAi-3 c t K«fc <9 . 

[0045] *mmtem<o^m$mx^ 

[0 0 4 6] *i\ xT7 n-;*-* 4 1 CO^mif-^-tC 

SMM l 2w«H*»AS«^iiia-rsiBlcWi«y» 
4 1 2^e>*^nSf^fi^PRA^S(^CT^t:i-5 

->-*4 leottaft^tRAffijRitoBBffli^jt 



*4 lco^a3ff^^^co||^y>^<i:^S<3v>T©A^ 

[0 0 4 7] Ls&^bftiSfe, _LfS;tf^;*Jf 4 1 4I3\ f& 
§4 1 2 fcJt« UT— Jft^JK»*3&S*#V^fc«>. -^COM 

KTgstt«ltBffi4 1 0<05£ffSriiigi-SKAS«O*co 
Affl^dSKftLftv^flfiBEWlcfco-Ctt. ±EBB» 

30 sc^s^coco. ra«AS«;**s*ftt--6iaatWEfiF«cfc 

oTIt B8*^^^e>*:a?>fetL5KA^fi[*^oia 

[0 0 4 8] ^IT\ *jafi»tt-C«t, ^ 9 ^fclf^* 
MA 1 4CO?^^>fb(COV^rcO/S«iltL^^*aSUT 
J£coB£A^fi£*^Sfc*6t;:. ^§4 i 4(cBli"S^ 
spflFe-T* /uSr** U coiS^e^vW-S^v A^« 

40 »C>ix5»A^«*S:^^^»[3(ffiGAi: LTKSJt 

So 

[0 0 4 9] £i\ 1 2 0i&«S: TThtj ftZ 

t m2&m) . !»»4 1 2^^7^14 1 4(C^jg 
^tiS^ffl:^. 1^14 1 4d^KA2SSt-ea*tt-5f» 
*i:<^«i[9-&v^^ x ^/7^®4 1 4co^gT gs co^t: 
id (Tgs) /d ttt. JL^TC05^ (1 8) coip^-T 

[0 0 5 0] 
[»12] 



(7) 
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11 



12 



•08) 



r^T\ taj , tbj , rcj , tdj ft. ikgfu l 

<b!f*^4 1 2 i<ORB0>j|&ea*. #^:*jf 4 1 4 £55 A 
[0 0 5 1 ] £fc. ^7^1414 

t. HUM l 2S.t5KA^mir^ra^fK^)SS^Si<^ 
5fc<£>. (19), (20) {c^-TJ: 51-, a (l 

8) <£>£iz2, m-htf^^m* i 4Mtatd (t 8 
s) /d ti* roj jc^d. &tz, T^ygsiGA^: 

[0 0 5 2] 
»13] ■ 

o=^ (r A/ -r s ,)-(c+o.^.(r A ,-r^) -a*) 

04 = mt -"(20) 



lot, ±ffi#^: (i 9) 
i) tmbtiZo. 

[0 0 5 3] 

mi 4] 



(2 0) ^feftO* (2 



04 



"(21) 



Mtc. (21) Srfflvvoe (is) Srf^ts 

2) ^e>H5 0 ^ 

[0 0 5 4] 
[Rl 5] 



dt 



= {^^ ^)\^-4GA) -(22) 



:^ (22) «oa£iaii, KftPflPB* rAtj £ 

(2 3) cOi9tC^i-r t36S-C#5o 

[0 0 5 5] 

mi e] 

= s — — *** (23) 

Sot, ±I5*S: (2 2) , (2 3) frb&LTom&iZ 
10 (2 4) ft. ^<£>5£ (2 4) 

ftT*£>9, £fc. S3Hrt<o r j j te^HMS. Tj-ij 

[0 0 5 6] 
1*17] 



a 

a7 



U - i) I 
» J 



-(24) 



r<D5t (2 4) ^P>M«b^*J: 9 SI*jSTO^P3/ 
20 h/U#iiiB:i:mt(j)rt\ SftSWv^/7 P ailGA 

(j)£, IKKIHRKA ttto^y^fltlMlCGAO-DfclcS 

[0 0 5 7] Wi/^ny b/U|a«TA^(D«BBaHE^ 

©*S.(2 5) , (2 6) HJ;oti^W 0 
[0 0 58] 
30 [ft 1 8 ] 



-(25) 



2k 



2 K 



a$#) oaawwiu rpaj i**»e {-mm -cs> 

5* TAtJ ft. 7D> 7 h;U#2 4Wf|PiiT* 

[0 0 5 9] h/UKSTAdS— 
StL5^*51<E«P*:«*i-Si:. * (2 5) (2 
7) (DJ:5i^^ 



[0 0 6 0] 
[*1 9 1 



•••(26) 



fff f7 a - * • A, ■ - 



TmtTAj % TPrnTAJ fit, 3£*iHE»^*5*t* 
SO -amtTAXtKTiKflJK^EEPmTAH. ±8tt«j«|5# 2 7 



(8) 
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13 . 

[00611 (27) £/BV^<t. * (25) tit 

(2 8) ©JipKlSt-i^T'tS. 
[0 0 6 2] 
IK2 0] 

m t = m tTA ZTT^ K**) 



14 



A£fS^KA^^ftmc£j&S^L< ft 9 (mtTA=mc) ; 
£7c. TffifflK««aTm^_h8KftiJ»S[?aThai:^b< 
ft5fc«> (Tm=Tha) % ±15* (14) 11^^ (2 
9) ^)J;3l^tr<i:^"Ct^ 
[0 0 6 3] 

HR2 ll 



m tTA 



0 " P mTA - *) 



•(29) 



tit, (2 9) £ffl<^5<t. ±1E5S (2 8) We 

(3 0) o«t 5^*"T-^d5"X?*So 
[0 0 6 4] 
[£fc2 2] 



r mTA 



■(30) 



w<D5t (3 0) ^bflf>^4J:5l^ ^oyh/u^iiia 
imtil S«F^T<z>±iffifl(K«?aThaRt^TiffiftiJ«ft 

t#OT8fift»J»^ffiPmOffiPmTA (TA, NE, VT) 

[0065] *jaK»HR-ett. ^(oxo \c&%\\(D^m^ 
(14) ~ (17) fcfc*<5v*TW*/£*»e>Brje 

i<D&mmmzm^xjfttbbtiz>xv v h/^mmm. 

[0 0 6 6] *IIJ6^ffi^»fi«cj:s^iJ«! ! 3ao 

3&T>*lKl4te. 1BrtKASaSmcS:^9J-f-5B8^>#JIIS 
m^^i-y d—^-v— h-e&>5o ilo^s 



*HE*«3 0^J:9gf^Bia^»i&Hffl («x.tf 

r 8msec. J fe) "TfH^^tl^o 

[0 0 6 7] ^(D&mizmisXlL KA2g«lftS:. 

T-£>5 0 Sot, ::m m^^bHuIe]^^*#t 
f-^ Lfcfl#Has&iS L^«Offlrt»A§gft*mc36STaiJ 

*acTl^m*^^sftSfco1a^rt®A2g«a:mc^c:s<5v^ 
[oo68] 3SNF^^e>^-aflfflratfwd«(DT«ftfis 

KSUEPnu TSfflHBWW&Tiiu ^T>*-€:ttbOJtPm/Tm 

tmrnx* y h/is&mmmmTzmnituzm^^xwftz 
iiia**^fimtnc:s<5v>r*asKs^»jffliajtfwd» 

0>T8fc«!IR»BEPiiu T*«»«flTii, ^^tt^coit 
20 Pm/Tm^r^tb^V TPmltfwdJ % TTmltfwdJ x 

rPmltfwd/TmltfwdJ £ s ^T^^tUa^^T^T 
«E<W»«fiEPin. TSEfi(K««aTiD. RTf*:ti(b<DttPm 
/Tm^, rpml(i)j . TTrnKDj . rp m l/Tml 
(i)j ^-ttt-t*tb*|Hi-S« 

[00 6 9] rOT*«»«JEPmltfwdXt5T*«l»« 
^TmltfwdO^mtC^LTfi. _hf£#5fc (15) ^ (l 
7) XO« (3 0) {cS<5<»a[S:Sf^lH]«»t»^5g 

[00 7 0] ip*>, ll4lOTt J: 51-. Si\ ro^i 
30 0H3Eia*Sr*-r*O^«[i («JSMI: roj ) j&W 

>^y^>h^^ (^7^202) . m^mmmtf 
wd&mjzmm a 1 1 «±-c*> s a^sa^ww £ (* ^ 

'7^2 04) 0 ::t\ ilEBf^ttWIAtltt, ^ 
>^#IiA^^^^U< («XLfi rsmsec. j ) SSJ^Stb 

tv^ 0 tLT, ^8HJ(lRltfwd36S^cDBfftl»BIAtieUi 
Xfo'b b.^\Wi t£ ttfz.Woi~i (^T 7 7^2 0 4 : YE S) . 

(is), (i6) <DMmmm& t&mmmm 
Mitm^<mm^ti^> (^^^206) 0 

[00 7 1] JSfetC. * (3 0) (^So*V>T^n 7 

40 iiig*ie3e«mti(i)3ftsimj*ix, M(-^com. 

(15) - (17) tdS-^V^T. T*«»»EPml 
(i). TaE«K«aTml (i) . &VM:tlh<Dtt Pml (i)/ 
Tml(i)^n^^fcfcS^nS (^7^/2 0 8) 0 ^ 

lt, ma^»JWiBtfwd^e>3fs«FH]Ati*«»$ 

*U toWHI (tfwd-Ati) j&sSffc^^aiflBMtfwdi: 
[0 0 7 2] J^tC. rofr/c^^iij^^tfwd^ roj 

50 2 1 2: NO) ^y"yzf2 0 2 £lfe<D%}M&ft8£ 



(9) 
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t f vdAmj&$ m a 1 1 *as T'frzt mm zinzt ( * ^ y 

-72 0 4 : NO). , BSffcffflPgA t tf^jHJHIHtfwd £ ^ L 

[0 0 7 3] -t L~C X ^77^2 0 2 — 21 2 <D^^^ 
^ll9ig^nsri^J:0, ^SSMBMtfwd^ roj J£*TK 
(^^2 1 2 : YES) . TStflSKSUE Pml 
(iK Tiffiffl»ftiB.Tml(i), ^t^lt^n^OPml/Tml 
(i)riv ^9J«BPfltfwd«^T8ttffllKft£EPmltfwd. TSc 
{Me&^tfSTmltfwd. &tW;Pmltfwd/Tmltfwd<t IT^ 

[0 0 7 4] #J;tfcf. ^^(D-m&mMfr rsmse 

c.J 9t<D^ s f^'fl 0 O^^^riiCT^SiJ^ 

Ppltfwd^S r3 0msec.j KK^£*L*:*g\ ftl&P^RSA 

t^±ia^^^/7 P 2 0 4, 2 0 6, 2 0 7^1^IC 
r 8 msec. J — » T 8 msec. J — > T 8 msec. J — ► T 6 mse 
c.J £V^o?d§^TilKK3££*U 777^2 0 2^ 
77^2 1 2^T*co^^^ T4 1hIj *95S§*x5r £ 

^tr, #*>*fi£ias r 4 j -cfca^g^TSE 

WKftEEPml (4) . TSEWK^taTml (4) . &t^E: 
tl^OitPml/Tml (4) *K ^Bfr&a>k^»JMHItfwd 
«OTJSS<B!JRSUEPraltfwd, TSEft'JK«iS.Tmltfwd. 
rJttPmltfwd/Tmltfwd,*: LT^*L^:ft,R3£$jft,3 r. £ 

[0075] ^rom, a* roj |:!/tyF 

L<R£*ix*:« (^^s/^2 16), AattBi 3 tc^ 

[0 0 7 6] ^f^/SOOm x77n-^^ 4 

(j) (but* tga(o)j tmttz) &&m£tiz>o & 

(«T* tga(-i)j tm&'tZ) tl-S^t. 5te<05£ 

(2 4) SrfflV>rS,^FJS"C^^oy h^#ii®*»J£fii[ 
mt2(o)^aj£H& (^7^^400)0 

[0 0 7 7] tLt, Z<DX*<% h/lsfrmm-MLfflfemm 
t2(0){cg-5^ ft<D«-5t (15) - (17) SrJfll^T 
KIRIHHfA t^TMMEPni(i), T««»«i&T 

m(i), X^^^<OJtPm/Tm(i)^tt^nifUi*^ 
6 (^7^5 0 0) 0 ft, rcOX^^/^S 0 OfC&l^ 
T#ffl£ft*Tffi{M»*UEPm(iK TSKffilRSHaTin 

(i)> S:t>^^XbOitPm/Tm(i)(7)9^ x 4*lE]tf>»i&* 
<>< Z^yX'&mZtlZm&ttl^tl rpm2(l)j % TTm 
2(1) J . rPm2/Tm2(l)j t , Su [3 <£>»]&*< ^ 

■e«asiLfc<l[Sr*tb-eix rp m 2(o)j . rrm2(o)j s 

rPm2/Tm2(0)j £^iEi-£ 0 

[0 0 7 8] ftd, £JlTa>#S; (31), (32), 
(3 3) «iS<5V>xm^^?>^»J«BratfwdS<OT8£«l 



16 



KfiCJEPmtfwd, TSEffilJK^Tmtfwd. &t>*-£*l, b<D&: 
Pmtfwd/Tmtfwd^gttB^tL-5^ £^1-. ±IS^ (1 
4) 35>e>»t>ttS* (3 4) (dS-^V^TTWMBratfwd^ 
*Ci3ttS«rt«A^fiCftmctfwddS|mi$ix5 (^7 7 
7^6 0 0) 0 

[00 7 9] 

It23] 

Pmtfrd = ^Iz/W - Pml <°> + P»»2 «» ' ' '(31) 
1U tntfwd 1 m\tfwd 'ml I m2 



40 



l mtfwd 



•'(33) 



=4^- • r- (« • ^ - *) • • '(34) 

r mtfwd 



ZL 0 UTfSrtKAS«*mctfwd^»ffi*^y5:« % ftlal 

R^^tL5 ^77^7 0 0) o 
[0080] H5tt, ^ 9 LtlO»I<Og^if|:/Tt 

[00 8 1] I^il^t«t5l^ 2*lHlB£JU»«>«ii^>f 

(^jt 2, t 3) ic&tfztnmm. bp*>t*« 

Pi^m^Pml(O), Pm2(0). TMW&^tELTml (0) , Tm2 
(0) x r^U^COitPmlXTmKO), Pm2/Tm2(0) > ^rj^ 

^xmmm^tik(om (Pmi (d ^vM^m^p^ 

■e^ffi (GA(0)) ^ LTStU^tLfcM^rffiV^^n^ 

[0 0 8 2] 
[§24] 

Pml (0) P^l (1) 
r«l(0) «- r«,l(D 

7-^2 co)- r W 2 0) 



(0) 



^1 



^-(0)4-^0) 

y m2 7 m2 

S^cofiiiZ.^^ W*Jt 1) \zj$rtZ>ti)M 

[0 0 8 3] 
M2 5] 



50 



(10) 
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Pmli Q >+- Pa 
Pml (0) - Pa 
T m \ (0) «" ^ 
r OT 2 (0) - T/rp 

i in2 y Aa 

[0 0 8 4] £x±.mw^tz£o\c *mmmmxn. t 
^mm^m^m^^x^^^mm u 

[0 0 8 5] ( 1 ) S£oT. T^JP^^Tm^^^ 

[0 0 8 6] (2) £fc. ^ 9 Lfc^A^m<7)^#^ 
®ft^IiJ^^^^/i-=¥i^#IiJ^3feo^; (12), 

(13) ^^^n^x^^mm^zx^xm^. :m 

<bl§t»5£ (12), (13) (Cg^<^^^[Hl^ 
t f wdf£ CO TSEftd B&fUE P mt f wd& t^T^J T mt f wd 
%&rnctfwd(D^m\Z&m?J:Tffimv&^I£ Pmtfwd& TJT 

[0 0 8 7] (3) MtC v ±fS^^^^O#J^ 

nxmn-tzmi^ &®(Dm&?4 ^^^xm^^btiz 

^H$;£T*<?>T^{R'J@mJEPm 

5 e r<Dfc£>, T^<JUi35#2 8^0/E^^^^tU-r 
ni(0)S.U?T?5fE{ail?^^Tm(0) ^±fa^ig^^ii CXI 

^^is r t ^ctp i^fts. 

[0 0 8 8] (4) ^^T'OTfflEllileS^Tm 
£5 etptcftSo 



30 



40 



18 



[0 0 8 9] ^SSff^r-fi, * * 

n-^-^4 1 C0^tHm^S^V>T»tt5$ttS^o y 
[0 0 9 0] r~T\ r.tbf>^o 3/ h/u#iiiga«l^ffi 

b^3TM^^U£Pml(i), Pm2(0), Pm2 (D&I/T 
10 gSWJK^jaTmKi), Tm2(0), Tm2(l)fi v HBSOlKft 

[0 0 9 1] gP"b, ^05/ h/u#iliSS«l^fflmtl$rffi 
V ^T*£> bttS TSKfflKftJE P ml (i) . TflMKttiBTiii 
Hi). &^tte><OJtPml/Tml(i)tcBIUr(i. 4* 

i8*ffl*fimt2S:ffiv>-C*aC)P>*L-5TMffilK«EEPni2 
(0), Pm2(l) v T«£fflIK«iaTm2(0), Tm2(lK IkXf^ 
20 *L£>(£>JfcPm2/Tm2(0), Pm2/Tm2 (l) icm bTf3\ ^ 

[0 0 9 2] Z.CDfz.&b^ $t<n%r3Z (31) — (33) £ 
fflV^T*»W*^fflSBWtfwdS<7>TStfllJK^JBEPmtfwd 
^t^T^K^MTmtfwd^attH-r^Cl^jC^D. iB« 



[0 0 9 3] 
[»2 6] 



PmtXjwd = PmiW 
P mt\fwd _ P m \ 



l mt\fwd 



(0) 



[0 0 9 4] 
[t27] 



P mtfwd 
T mtfwd 



Pm2 <°> 
^2 



(0) 



'm2 



50 



<h ft 6 n ^ ^ ^> . TSSffi9»«/E P mtf wd&t^T«£fiBK« 
^Tmtfwd^^ n v h/U#iiiS*»J^ffimt 2 (7)^ldg<5 

[0 0 9 5] (5) ^oT. *jafi««^J:ntf, 3£« 

wd«OTffifi!lKftJEPmtfwdXr/Tffifil]Kft?aTintfwdSr 

ctfwd&J: 9 jE5t^^W-rS-i:^-Ct6«t 5icft5 0 
[0 0 9 6] (6) JEtC, *Hlfi?gffi-ett, 



(11) 
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20 



* 4 1 (0^7^14 1 4 iCitei^f/U^i 
aij^fimt2^*^5 J: 9t-LT^6*:#>. ^14 14 

[0 0 9 7] 2 tomtiSMffl UT. 2 O 

[oo9 8] *nmmmx-i*. «A^ft**:«ia-rs-ir 

ttt. (H6^i-J:9tC s -fr— v?^ 2 2{dR0f**te> 

[0099] *xibeiRoass(c ± s^aggy®^ 



(#J;t{3 TSmsec.J S) T^^T^tto e ft. H 

ari^t^fy^i o oR^fy^2 o om 3fe 

cog] 3 RZ/m 4 (d^-r ^ 0.0^x^/2 0 

o—2i6 tmm<otimi>mn£tiz>tz#>, tco^m^ 

[0 10 0] rtl^^ryyiOO, 2 0 0(7)#M^ 

io (j^t. r.Pm3(o)j £*iei-3) **»ui*ns. & 

a*4x^T«K«fl®SUEPm(i-l) (J£*T. rp m 3(-i)j <b 
SIM'S) 5t« (15)^f>#^5W (3 

5) ^Kg^vT. ^jftd^MIIRnniA t ittreo (gp 

(i-1) (J£AT. r^u 5, h/U#iiia*»JSffimt3(-l)j 
t^0) &&m&tlZ> (^7^4 5 0) . 
[0101] 
20 [»2 8] 



m /3 (-l) = 



mt3(-l) 



-*>} 



--(35) 



(i6), (17) i^f'bnsftoa (3 

6) , (3 7) t^^S^T, W»jS-C©TSE«ft««B. 
Tm(O) (J^T. rTm3(0)j tmm-fZ>) % afcUHCTSt 38 



^W^m/EPm3(0) <tT««»fti&Tiii3(0) k <Dit Pm3/T 

[0 10 2] 
2 9 ] 



f * ( m,3( " l) * ^tU * (a ' ^ 3< "° ~ b) ) " * (36) 



7 m3 



'(37) 



ft. -blE^T-y^ 5 0R^7^^5 5 OtCjo^T. 
5£ (3 5) Xt5 (3 6) |cS<3V^TT8Kfl!lK««aTm3 
(0)^£rJ?tM-S#<3\ ±*«R*kThaf±1ttH©Ma i 

sm^lcf*. ±ffi{BI»«jfiTha«*S^.S"C-(Offi. 

a, b»«n«Jk«O«0BIsMEfiCNE (= r 0 j ) RXf 

[0 10 3] Rtd, jfcp** (3 1), (3 2) ,f£:fc5^ 
Tx rp m 2(0)j , rp m 2XTm2(0)j Sr-ttL-ftt rp m 3 
(0)j s rPm3/Tm3(0)j KBJfeLfcsti: . «-5t (3 

3) , (3 4) t&m^x* m$&frt>^-mmmtfwd& 

(OTSEfilK^ffiPmtfwd. TSSffllKftaTmtfwd. 
tt^OWlPmtfwdyTmtfwd^^tfci^tL^i: <h *bfd N 50 



aratfwd«(c*5JtsiarrtsiA^*mctfwd^a[tB$^ 

5 (^77/6 5 0) o 
[0104] r 5 Url8rt5RAS«*mctfwd3ftS*tb$ 

ffi^^^t-etLias^ixs (^7^7 5 0) o msfs. 

[0 10 5] |^H!^n5l:; 2IeJS^I^OfiJ^^^ 

»at 2, t3) \a$tfz>$}mm. w-^Tmm 

eftMJEPml(O), Pm3(-1). Tatffll»«ig.Tml(0) f Tm3 
(-1). ^T>*{crtbbOM:Pml/Tml(0), Pm3/Tm3(- 

ffi (Pml (1) ^V>H^^X-^ffii LT*^^tL 
fc« (Pm3(0), Tm3(0)^) ^rffi V ^ tth^ftl^TCO X 

[0106] 



(12) 

21 

[§C3 0] 

Pm\ (0) «" (0 

-Sr. *3O0>9J&*-f (B^Jt 1) {C^ttS^ 

[0 10 7] 

Ii£3 l] 

P«! (0) <" 
PmSC-D <- P a 
T m i (0) - 7V» 

T m s «- 7*^ 20 
Y m3 

m 1 Lfc (1) , (2) . (4) 

m i o^^^teSc ufc (3) ^fflMco^t 

[0 10 9] Mt-. P^iEir>^4 6 O&fcBfl ^(clo* 
V^T^^^tL^Tat^JP^ii Pm3(0), Pm3 (1) . T 
SWJH^&Tm3(0), Tm3 (l) , &tWx ib<£> it Pm3 40 
/Tm3(0), Pm3/Tm3 (l) t± x ^i&EB^Cljo^T 

a-r^M^^fc^o ffiot, ^mmmmA^x^xhrni 
<DnmmmAc&\,^xmwiL,fc (5) if^^of^^m^ 

[oiio] ^mxistz&mmmmit, £it<dx?\c 
mtizz^ i,xmMi-z> ^thx% z> 0 
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m-5^xm=f-m*?r o x 9iattJ;v\ 
[0111] - ±is^^^9ir'fi. pjj^# 2 e <ds<ju 
^y&m^zmm&mttftmmm 1 o 

[0 112] • £fc. ^Jx^^ilK2 O^A^O 
«ft«r#J»+S«*. ^ix^^ 2 9 fCf§£i-3 * y 

m^^^-^commm <Dm$kti,xwt7£~tz> 
[oii3] • ±m&mmmm<D^®mmzio\,^xte. 

Pi^^Tm^t^T«W^iEPra(D^S:^i-5|l^ (1 
4) £r/8^5 J: 9{cL7c^ ^^ftmc$r*:i6^^ 

[oii4] - ±^mmmmx^ ^A^m^^^-e 

br#^tL^_hta#^: (8) . (9) zmmt 
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